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Effect of Radiative Heat Transfer on Cosmic-Ray 
Transport in a Rotating Cloudy Interstellar Medium 
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Abstract — Radiative heat transfer on cosmic rays (cr) in a 
rotating cloudy interstellar medium is modelled by imposing a 
time dependent perturbation on the cosmic transport in the 
interstellar medium containing randomly distributed giant 
diffused molecular clouds the temperature involved are assumed 
to be large so that radiative heat transfer is significant. This 
renders the problem very nonlinear even on the assumptions of a 
differential approximation for the radiative flux in a adiabatic 
concentration and zero activated energy. When the perturbation 
is small, the transient flow is tackled by the laplace transform 
technique, giving solution, for steady state spectrum of cosmic 
ray in the galaxy temperature up and down stream 
concentration, also velocity analytical and numerical solutions 
were obtained for temperature. Velocity of the down and 
upstream concentration in cases where radiative heat transfer 
through the rotating medium are taken into account. The 
incorporation of radiative heat transfer is of particular 
relevance to cosmic ray transport whose temperature is usually 
very high which radiate a lot of heat through the media which 
surround them. 


Index Terms — Cosmic Ray, Radiative Heat, laplace 
transform, interstellar medium. 

I. INTRODUCTION 

Modeling of cosmic rays(cr)transport in the cloudy 
interstellar medium containing randomly distributed giant 
molecular clouds have been considered in the works of(l) 
cowsik and wilson(1973) dogiel et-al (1987)monfill et-al 
(1985),(2)osborne et-al (1987) and ptuskin et-al( 1990). the 
latter two papers the theory of diffusion in the cloudy medium 
which takes into account the cloud’s finite transparency for 
diffusion particles. Due to the high temperature involved its 
application in astrophysics, meteorology cannot be 
overemphasised .In all these investigations the problem of 
radiative heat transfer has been ignored unfortunately high 
temperature phenomena abound in cosmic ray transport 
which radiate a lot of heat through the media which surround 
them. A primary difficulty in thermal radiative heat transfer 
studies stems from the fact that the radiative flux is governed 
by an integral expression and one has to handle a nonlinear 
integro-differential equation. However, under fairly broad 
realistic assumptions the integral expression may be replaced 
by a differential approximation for radiation. Thus in one 

space coordinates z, the radiative flux Q satisfies the 

nonlinear differential equation cheng (1964) similarly 
arrienus will be linear zed. 
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Where T is the temperature of the CR C is the 
Stefan-Boltzmann constant and CC is the absorption 
coefficient which will be assumed constant in the modelling is 
the chemical reaction. Taking into account the cloud’s finite 
transparency for diffusing particles CC and equation l.laand 
1.1b is approximated by 
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dZ' 

— Fk — 
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in which subscript oo will be used denote condition in the un 
disturbed cloudy interstellar medium. In the subsequent 
analysis, the mathematical formulation of the problem for 
finite transparent medium is presented in section 2 for the case 
in which the temperature of the medium suffers time 

t' perturbation when this perturbation is small, and the 
non-linear system tackled by asymptotic approximation so 
that the first order transient problem is solved by the laplace 
transform technique the results are obtained in section 3. In 
section 4 a discussion of the results of the previous section 

given. 

II. NOMENCLATURE 

i 

q _ = radiative heat flux 
q = complex velocity 
k p = Boltzmann constant 

— Stefan-Boltzmann constant. 

,2 

H o = constant transverse magnetic field 
= thermo diffusion constant 

1 2 

k r = constant associated with chemical reaction in the 
Arrehnius term 
P = Prandtl number 

r 

R = Radiation parameter 
D f = Diffusion parameter 

i - V— T complex number 
E = rotation parameter 
M = magnetic parameter 

G c = free convection parameter due to temperature 
D n r = mass diffusivity 
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S t = Soret parameter. 

S c = Schmidt number 

(u , v ) = dimensional velocity component 

(x , Y , Z ) = dimensional Cartesian co-ordinate 

k = thermal conductivity 
g = gravitational acceleration 

c p = specific heat at constant pressure 

T = dimensional temperature 
C = dimensional concentration 

= reservoir temperature 

= reservoir concentration 

T aj = constant plate temperature 

c co = constant plate concentration 

T m = mean temperature 

G c = free convection parameter due to concentration. 

}.i = is the permeability 
( J = absorption coefficient 

,u = coefficient of volume expansion due temperature. 
d c = electrical conductivity 

Y = coefficient of volume expansion due concentration. 
E= is small parameter 
IT = kinematic viscosity, 
fi = angular velocity 
pm ^reservoir density. 


III. Mathematical Formulations. 


We consider the cosmic Ray transport (CR)in the interstellar 
medium containing randomly distributed giant molecular 

clouds with characteristic velocity U 0 and rotates about the 

X-axis with angular velocity Cl .This is maintained at 
temperature and 

.concentration. Td + £)/( t') And C K (1 + s)f ( t') In 


which, T k » 1 , / ( f ') ,C K » 1 


/( t r ) is an arbitrary function of time and £ is a parameter is 


small. Following the arguments in Opara (1990) and 
employing equation 2, 2b. The governing equations for a Q 
medium are; 


du 

~df 


2 f 


- 2Qv' = v 
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+ g/?(T-T \) + yS(c-c t ) 
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Where U=,U Q ,li =0 ,T=T k [1 -\-sf (t ' )] ,and 
C=C fe [H-£/(t f ] on Z'=0,u'=0 
l/=0,l/=0,T=r o3 ,C=C i;b onZ' -s 00 

[7] 

where (U',V',0) is the velocity component, k r is the thermal 


conductivity, g is the gravitational acceleration,^ is the 
electrical conductivity, ju, is the permeability, v is the 
kinematic viscosity, a, fj are the coefficients of volumetric 

expansion for concentration and temperature, C ^ is the 
specific heat capacity at constant pressure, C is the up and 


down stream concentration and D r is the diffusion 

coefficient. For simplicity we shall take f(t * ) = 3~C ^ )the 

Heaviside step functions. 

Introducing the following non-dimensional quantity. 


(0. 8 k )= ( C f C k ) = (C.C k )/C QD ,E=nv QD /V.l, 

t =t'ulfv ,{U r V)=(U T V f )/U 0 
M 2 =a c (ji z K 2 /p aj ,G r =g(fi , u ^ , 

G c fiL p J C co ■ 


[ 8 ] 


Acer, 
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Equation( 3) to (7) can be written as 


+ 2iEq — — M z q + G r {6 — 1) + 

G c (C - 1) 
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[9] 

[ 10 ] 

[ 11 ] 


such that 


Z = 0; q = 1 , 6 = 9 X [1 + sK (*■)] , 

C= CjfC 1 + ?K(£)\,Z m: q = 0; 8 = 1, C= 
1 [ 12 ] 


and Cf — U + ffi — V — 1 ,E is the rotational 

parameter,, “ is the magnetic parameter 

G c and G r are the free convection parameters for 
concentration and temperature respectively. 

is the total radiation parameter. The mathematical 
statement of the problem embodies the solution of equation 
(9), (10) and (11) subject equation (12). 


IV. METHOD OF SOLUTION 



The problem posed in equation (9), (10) and (11) very 
non-linear and generally will involve a step by step numerical 
integration by (say) the explicit finite difference scheme, if 

however £ is small analytical solution could be possible by 
adopting regular perturbation. A great physical 
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insight could be shed on the whole of the CR transport by such 
a perturbation scheme and this is the problem we shall pursue 
in this investigation. We write: 

8 = 8 m {Z) + £e ( 1 ) (Z,t)+~ 

c= c co: '(£) + £C a H_z,t) + - [i3] 

q = q' : '°'(Z) + £q ( - rj (Z, t) HI — 

Substituting equation (13) into equation (9)— (ii ) we have 
the sequence of approximations equations. 

2 

4Eq = " ■M‘V G) + G r (8 m - 1] + 


dZ 


G c (C C0) - 1) 


- 4 %(0 Co) - 1 ) 
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d 2 C° 
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[15] 

[16] 


dZ 

z = o -8^ = e K , c m = c K 

,q w = l,Z-> go; 8 m - 1; q c ° } - 0,C TO = 1 

[17] 

and 

+ 2 iEq 1 = - M 2 q m + G r 8 m + 


dt 


dZ 


G c C fl) 

[18] 

391,1 = 329 “ - 4S s .9 m ®“ ) 

dt ez a J 

[19] 

ac ci] d 2 c m m rn 

dt dZ * J 

[20] 


Z = 0:8 = CjfJ f? C1) = 0:Z -> 

oo; 0C 1 ) = 0: C'* 1 = 0 

= 0:* = 0: = 0 = C= q (1) :Z > 0 

[ 21 ] 

In equation (21) we assume 

£ > 

0 .Thus by asymptotic approximation ., 

the problem now splits into a steady flow on which is 
superimposed a first order transient component, the solution 

for 0WandC TO are given. 


z = VC-H-) C 




ds 


° ^ fs s -5s+4)2 
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[22]a 
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and 


q KlJ — exp [— (M 2 + 2'iE') *Z\ — 


V (^-J &o dz 


231 


smh[ 


M 2 + 2iE) i{z(e {o \ C -2(5}] 


Cy-l)ds 


(S 5 — 5s+4)2 


[23] 


To solve equation (18), (19) and (20) and also consider the 
boundary conditions in equation (21) we take the Laplace 
transform with respect to time. Denoting the transformed 

variable by f and placing a tilde over the transformed 
function, the equation satisfied by <7 8 < '- 1) andC < ' 1} * 

d2qW -{M z + 2lE')q {1: = -[G r 0 (1) + 


are: 


dZ 


2 


G c C 


(!)] 


.2 


dz' 
d s cW 

dz 2 


- (43l T 8 m + £> 0 & = 0 


- + 



= 0 


[24] 

[25] 

[26] 


with boundary conditions 
Z= 0,8& = 8 k ,C^ = C k and 

q m = -[G r dW+G c cW] 

z~* os, — c C 1 .' =0 = q 
The solution for 8 and C ^ are 


(1) 


[27] 


0« = 


exp[-{4%+|)sj] 

? 


[28] 


and 


= 


exp[-(;4%+f)2,Z'] 


[29] 


Hence, taking the inverse Laplace transform 


1 1 
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i 
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Next we consider when 8 k , C k are arbitrary and approximate 
8 (0 \CW by 

8 m * = (8$ - i)e~ 2fiz + 1 [32] 

_ ^3 _ yjg-ZYZ _|_ i [33] 

By virtue of equation (32) and(33) the solution for and 
now reduces to 
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[35] 


where 


r) = 3l T {B k -i) 

and 

J n (x) andl n (x) are the b Bess el and Modified Vessel 

function of the first kind respectively. Equation (34) and (35) 
have simple pole at c = 0 and another branch at 

f = 4% and R - respectively .These Equations could 
then be inverted in the (39) and Bromwich contour with a 
suitable branch cut and we obtain 
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Because of the complexity of the above equations (36) and 
(37) we find it expedient to consider limiting values with 


/]jW 




~ i (ff'H n 

(2 Tirt) 2 


-> 00 


[38] 


We can show that 


l^lr 


Cl] 


-JC 


ft* 


exp [- 


j3Zf4J! r +£)5 

f 


] 




y3(fl r +^)z 

exp [- ] 


We the condition that f — > CO is the form of equations (28) 
and (28) we can now write. 


^-{e 2 P R T z ei'f c[-^~t — (4j r ? T t)z] 
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and 


When R t is large and of order 0(1) then 

» w fr (PA 

erfC — I And 


t 
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Also for -fl- — * 0 ; fol lowing Watson ( 1981) 

l n (V) = J 0 (x) - nK 0 (x) 

where Kq (xj is now the modified Bessel function of the 
second of order zero. Therefore 


e<Xi 

ft. 


i ^-Uoive-^) .j+ 

J D W < *_ 
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erfc(4R T t')2 ]j a s t -> 0 

fll r 

The solution for ■ ■ Ul equation (24) can now be solved 


putting 

S 2 = $ + M 2 + 2iE, 

We get 

q w = |G//e^ 

where 


0 (x) rfz 


[47] 
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LT -1 [e * J ]= e ~i M 2+liE ^iT-\ e C ] 

- e -i M2+2i ^— Lj g-{3-3)2/ t [48 ] 

and by convolution theorem we can deduce that. 

(,C1) = _L f‘dz f t J_ e -(«- , +2 i OT e -(2-r)/« x 


2jt 




'0 ± 
T2 


[(3 '■ 1 •' (t — z'JdT + C 1 -' (t — T)rfi] 


The solution is now complete. 


[49] 


Fig 1 : plot of variation of concentration and Temperature with 
Z for R is 1.0 
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Fig 2: plot of variation of concentration and Temperature with 
Z for R is 0.5 



Fig 3: plot of variation of concentration and Temperature with 


Z for R is 0.25 

V. Results and Discussions. 

The problem of the effect of radiative heat transfer on cosmic 
rays (CR) in rotating cloudy interstellar medium has been 
solved making fairly realistic assumptions, for a small 


time-dependent perturbation of the temperature and adiabatic 
concentration, the non-linear problem is tackled by 
asymptotic approximation, giving solutions for steady state 
flow on which a first order transient component is 
superimposed. 

Equations (22) and (23) give solutions for the steady-state 

components of the temperature adiabatic 

concentration c<*> and velocity fields. Equations (22) 
have been evaluated by numerical integration using three 


values of Rf (0.5, 1.0 and 1.5 to show the dependence of 

8&\C m on Z. These are shown in figure 1, 2, and 3. It can 
be seen from these curves that their cube values say 

8 m2 andC (m approach unity asymptotically. 

Furthermore, when the parameters are y = _ 0 


and /? = 12, the variation of the adiabatic 

concentration^ ~ and temperature 8 * with Z for steady 

state flow are given by equations (32)and (33) These are 
almost identical to that deduced from equation (22a) and 
(22b) .They are illustrated by the red line dashed curves in 
figures 1,2 and 3. 

Also for the case where 8 ' L — 1 the transient component of 
the temperature exhibit a standing wave structure. The 
magnitude of the standing wave increases intensely when the 

solute concentration C (Q) = x 

.This means that 


for 8 1 — 1 — the intense radiative heat is large, 


that is when R j is large. There is contrary to a relatively low 

radiation when Rj is of order 0(1) in which the standing 
wave structure is no longer exhibited by both the temperature 
and concentration fields. These are shown in equations (41) 
-(43). However, equation (49) give a complete solution of the 
transient component of the velocity field. 
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